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Abstract 

Interference between the halves of a double-well trap containing a Bose-Einstein condensate is 
studied. It is found that when the atoms in the two wells are initially in the coherent state, the 
intensity exhibits collapses and revivals, but it does not for the initial Fock states. Whether the 
initial states are in the coherent states or in a Fock states, the fidelity time has nothing to do 
with collision. We point out that interference and its fidelity can be adjusted experimentally by 
properly preparing the number and initial states of the system. 
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I. INTRODUCTION 



Since the recent experimental realization of Bose-Einstein condensation in small atomic 
samples , there has been much theoretical interest focused on the physical prop- 

erties and nature of Bose-Einstein condensed systems such as coherent tunneling, the col- 
lapses and revivals in both the macroscopic wave function and the interference patterns 
ffflflflflHQQ. It is hoped that the study of those exposal sys te ms will 
give new insight into the physics of BEC. Since the current understanding of BEC is largely 
influenced by the concept of a macroscopic wave function, the study of this feature is of fore- 
most importance. The investigation of interference phenomena should be perfectly suited for 
this purpose. Another motivation for the study of interference properties is the envisioned 
development of a new source of atoms, based on BEC, with high flux and coherence. It is 
expected to stimulate atomic interference experiments. 

Recently, Javanainen et al. have theoretically studied atom number fluctuations 

between the halves of a double-well trap containing a Bose-Einstein condensate, in which 
the two-mode approximation is used, which assumes that only two one-particle states are 
involved. They have developed an analytical harmonic-oscillator-like model and verified 
numerically for both stationary fluctuations in the ground state of the system, and for the 
fluctuations resulting from splitting of a single trap by dynamically erecting a barrier in the 
middle. 

This paper is organized as follows. Sec|Tl] gives the solution of model. Sec. IIHI studies 
collapses and revivals of interference intensity. Sec. II VI investigates fidelity of interference. 
A conclusion is given in the last section. 



II. MODEL 



In a symmetric double-well potential, the ground state of a single particle is represented 
by an even wave function ip g that belongs equally to both wells of the potential. Provided 
that barrier between the halves of the potential is tall enough so that the tunneling rate 
between the potential wells is small, nearby lies an excited odd state ip e that likewise belongs 
to both halves of the double well. The starting point is that we take only two one-particle 
states ip g and ip e to be available to the N bosons. The reason for this choice of including 



2 



only the two lowest lying modes for the double-well potential is that the other modes are 
energetically inaccessible. 

We adopt the usual two-particle contact interaction U(r 1: r 2 ) = (4-irh 2 a/m)5(r 1 — r 2 ), 
where a is the s-wave scattering length and m the atomic mass. Given the restricted state 
space of precisely two one-particle states, the many-particle Hamiltonian is [13] 

H = i(ei + e 2 )(a^ai + a%a 2 ) + ^(ei - e 2 )(a] l "ai - a£a 2 ) 



+Kiia,i 2 al + K 22 a\ 2 (X^ + Ki 2 (a^ 2 a\ + a 2 2 a\ + Aafaia^ a 2 ), (1) 

where we set H — 1, and correspondingly use the terms energy and (angular) frequency 
interchangeably. In Eq. (1) a x and a 2 are the boson operators for the excited and ground 
wave functions. The constants e and K are the one- and two-particle matrix elements 

ei = J rf 3 r^ e (r)[-^-V 2 + V(r)]^ e (r), (2) 
*2 = / rf 3 r^(r)[-^-V 2 + \/(r)]^(r), (3) 
K U = — fdh\Mr)\\ (4) 



m 



= — 

m 



J d 3 r\^ 9 (r)\\ (5) 



K 12 = K 21 = ^-J rf 3 r|^(r)|V 9 (r)| 2 , (6) 

V(r) is the symmetric double- well binding potential. Without restricting the generality, we 
assume that the wave functions ^ e>g are real. To simplify the discussion, we set K 12 = K u = 
K 22 = K, and e\ = e 2 = e. 

In order to solve Eq. (1), we introduce the following transformations 

ai = -L(A ie ikt -tA 2 e- iKt ), (7) 
a 2 = -j=(A 1 e ikt + iA 2 e- iKt ), (8) 
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where [A i: Aj] = Sij. We have from Eq. (1) 

H = e(AfA 1 + A\A 2 ) + K[(AfA 1 + A\A 2 f 



-3(A+A 1 + A+A 2 ) - 3A+A 1 A$A 2 + {A+A l f + (A+A 2 ) 2 }, (9) 
If we define two bases as follows 

\n,m) = -^=At n A+ m \0,0), (10) 
Vn!m! 

\n,m>= ^^af n a+ m \0,0>, (11) 

We have 

H\n,m) = E nitn \n,m), (12) 

with 

E nt m = e(n + m) + K{2m 2 + 2n 2 — mn — 3n — 3m). (13) 
We now define two-mode coherent states as follows 

\a u a 2 >= D ol (Q!i)Da 2 (Q!2)|0,0 >, (14) 

\u!,u 2 >= D Al ( Ul )D A2 {u 2 )\0,0>, (15) 
where the displacement operators are defined by 

D ai ( a i) = exp[a*Oi - a^aj] (i = 1, 2), (16) 
D Ai (u t ) = exp[u*A l -u l At} (i = 1,2). (17) 

It is easy to see 

1 i 
K,a 2 >= l~y|( a i + «2)e" iX *, ~y|( a i ~ ^2)e lKt ). (18) 

Considering the arguments of Bose broken symmetry, we assume that two condensates 
are initially in the coherent state. So that the wavefunction of the system at time t can be 
written as 



oo 

\m >= e~ N ' 2 V - T =(u 1 e^ K T^2e iKt r^pHE n , m t}\n,m), (19) 



where 

«i = — f={ai + a 2 ), u 2 = -j=(ai- a 2 ), (20) 



N=\a 1 \ 2 +\a 2 \ 2 = \ Ul \ 2 + \u 2 \ 2 . (21) 



III. COLLAPSES AND REVIVALS OF INTERFERENCE INTENSITY 

For convenience, we rewrite Hamiltonian (9) as follows 

H={e- , $K)N l + (e - 3K)N 2 + 2KN 2 + 2KN 2 - KN r N 2 . (22) 

where N{ = AfAi(i = 1,2), 2KNf and 2KN$ stand for two-body hard-sphere collisions, 
—KN\N 2 describes the collision between the atoms of the two wells. 
The dissipation is included by considering the master equation 0, Q 

^ = -i[H, p]+J2 H( 2A iP A t - A l A iP ~ P A t A il ( 23 ) 

J=l,2 

where jj (j = 1, 2) denotes the dissipation or loss rate due to some relaxation processes such 
as the coupling of the atoms in the two wells with the environment. To solve Eq.(23), we 
can introduce a transformation c = exp(iHt)c exp(-iHt) , then Eq.(23) becomes 

% = E 7,(2^/51; - A]Afp - pAfAi). (24) 

j=l,2 

The master equation (Eq.(24)) can be solved exactly for any chosen initial state. In 
particular, when the atoms in the two wells are initially in the coherent state (cs) |cki, a 2 > 
or in a Fock state (Fs) \n,m >, the corresponding density matrices are given by, respectively 



and 



where 



p( Cfi )(t) = |ct ie - 71 ', a 2 e" 72 * >< cue" 71 *, a 2 e" 72 *|, (25) 

n m 

p( Fs \t) = ^^( e - 27l ') n ~'(l - e- 2rrxt ) l (e-*» t ) m - k (l - e - 272 ') fc 

1=0 k=0 

xC l n C^\n - l,m - k >< n - l,m - k\, (26) 
Ax\ai >= ai\ai >, A 2 \a 2 >= a 2 \a 2 >, C™ = m!/n!(m — n)\. (27) 



The Schrodinger picture field operator for the sum of the two modes is ip — (Ai + A 2 ) / \/2, 
where the spatial dependence has been suppressed Q, E|- The corresponding operator for 
the intensity of the atomic pattern is and its time-varying expression can be obtained 
by the trace operator I(t) = Tr[p(t)if) + ip]. When the atoms in the two wells are initially in 
the coherent state (cs) |ai,a2 >, one has 

= - ^ |a i | 2 exp(-27 i t) + |aia 2 | exp[-r(*)] cos (28) 
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3=1,2 



where 



r(t) = ( 7l + 72 )t + 2 ^ \ aj (t)\ 2 sin 2 ~Kt, (29) 

3=1,2 

(f>(t) = /3 + (— l) J '|«j(*)| 2 sin5iir*, (30) 

3=1,2 

|ct/(i)| = \otj\ exp(— jjt), a\a 2 = |ai« 2 | exp(— i/3), (31) 

where we have set «i = |ai| exp(i0 Ql ) and a 2 = \a 2 \ exp(z0 Q2 ), [3 = <p ai — <p a2 . Eq.(28) can 
be expanded as the form 

/M W = 2 E \ a i\ 2 e M~2ljt) + \aia 2 \ exp[-( 7 i + 72)*] 
3=1,2 



00 00 00 00 

x 

m=— 00 n=— 00 p=—oo l=—oo 



E E E 53 j m (i« 1 (or)i»(i« a (*)i 2 )J r p(i«i(*)i 2 )Ji(i« a (*)i 5 



x cos{/3+ [5m- 5p + 5n + 5l}Kt}, (32) 

where J p (x) and I m (x) stand for the Bessel and modified Bessel functions, respectively. 

It is clear that when the dissipations are neglected (7, = 0, j = 1,2), and we take the 
terms 

m—p + n + l = 0, (33) 
we obtain a nonzero time-averaged value of the intensity of the atomic pattern 

/(cs) = ^(l«i| 2 + l«a| 2 ) + |aia 2 | cos/? 



00 00 00 

x 

m=— 00 n=— 00 l=— 00 



E E E ^m(kl 2 )/n(MVmWKIVKM 2 ). (34) 
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FIG. 1: Diagram of the time evolution of I^ cs \t). The time is in unit of x = Kt. The result is 
shown for the case of 71 = 72 = 0, when the total number of atoms in the two wells is N = 41 with 
|qi| = 5 and |ck2 | = 4. Here, (3 = tt/6. 

Eqs.(32) and (34) show that the intensity exhibits the revivals and collapses. This phenom- 
ena also can be easily seen from Figure 1. 

On the other hand, when the atoms in the two wells are initially in a Fock state (Fs) 
\n, m >, one has 

1 n m 

j(Fs)^ = i^^( n + m _ fc _ / )( e -27 1 <)n-« (1 _ e -2 7l4 y (e -2 72 t )m -fc (1 _ e -2 72t) fc^ C fc ) (35) 
l=0 k=0 




FIG. 2: Diagram of the time evolution of I^ Fs \t). The time is in unit of x = 71 1 = 72^ we have 
set 71 =72. 

which shows that the intensity does not exhibit collapses and revivals (see Fig. 2). 



IV. FIDELITY OF INTERFERENCE BETWEEN THE ATOMS IN THE TWO 
WELLS 



The fidelity and its loss rate of interference between the atoms in the two wells may be 
characterized bvjl3| 

F=< ^ |p(t)|Vo>, (36) 

L = -<M^o> |t=o, (37) 

where |?/>o > is the initial state of the system, p{t) and dp/dt satisfy Eq.(24). 

We now turn to study the fidelity of interference between the atoms in the two wells. 
When the atoms in the two wells are initially in the coherent state |ai,«2 > ; the corre- 
sponding fidelity of interference is given by 

F (cs) = exp[-|ai| 2 (l - e-^f - |a 2 | 2 (l - e"^*) 2 ]. (38) 

Similarly, for the initial Fock state, one has 

F {Fs) = exp [-2n7i* - 2m l2 t] . (39) 

Diagrams of the time evolution of F^ cs \t) and F^ Fs \t) see Fig. 3. 

1 

0.8^ 
0.6- 

F(t) 

0.4- 
0.2^ 




FIG. 3: Diagrams of the time evolution of F^ cs \t) (solid line) and F^ Fs \t){dash. line). The time is 
in unit of x = 71 i = -For simplicity, we have set 71 =72. The total number of atoms in the 
two wells is N = 41 with |«i| = 5 and \a 2 \ = 4. The Fock state is supposed in |n, m >= |1, 2 >. 



8 



Because the large particle number of the two condensates implies small ^jTFid, we can 
set 1 — exp[— •jjTFid} — JjTFid, such that 

F^expHla^ + WW], (40) 
for short time. The resulting fidelity times are then 

= [|«i| 2 7? + l« 2 | 2 7 2 2 ]- V2 , t<™ = (2n 7l + 2m 72 )^, (41) 

which show that the fidelity time is not only related to the initial state of the system, but 
also to the dissipation parameters. 

Furthermore, we can get the fidelity loss rates: 

l(cs) = Qj l(Fs) = 2n ^ + 2m72; ^ 

which indicate that when the atoms in the two wells are initially in the coherent state, the 
fidelity loss rate of interference is zero, but for the initial Fock state, L^ Fs ' is related to the 
initial particle number of the system and the dissipation parameters, but not to the collision 
parameters. 



V. CONCLUSIONS 



We have studied interference between the halves of a double-well trap containing a BEC. 
It is found that when the atoms in the two wells are initially in the coherent state, the 
intensity exhibits collapses and revivals, but it does not for the initial Fock states. The 
interference intensity is affected by the collision and dissipation, but for the initial Fock 
state, it is only related to the dissipation. Whether the initial states are in the coherent 
states or in a Fock states, the fidelity time has nothing to do with collision. For the initial 
coherent states, the fidelity loss rate is zero, but for the initial Fock states, it is determined 
by the initial particle number of the system and dissipation. This shows that interference 
and its fidelity can be adjusted experimentally by properly preparing the number and initial 
states of the system. 

It is pointed out that the recent realization of a superfluid-Mott-insulator phase transition 
in a gas of ultra-cold atoms in an optical lattice is very similar to the state preparation 
assumed in this paper, we hope our results obtained above will be useful to study Mott 
insulator phase transition in the future. 
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